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FIELD OF THE INVENTION 

The present invention is directed to a method of adaptive digital 
filtering and more particularly to such a method utilizing a reduced rank 
technique. 

BACKGROUND OF THE INVENTION 

Adaptive digital filters are known for use in various signal 
processing applications including speech and radar processing, adaptive 
beamforming, echo cancellation and equalization. These filters estimate 
the optimal filter coefficients from observed data, i.e. data representing 
the received signal to be processed. 

Reduced rank filters estimate the filter coefficients with a 
relatively small amount of observed data. Various reduced rank 
techniques are known including Principal Components, Cross-Spectral 
and Partial Despreading methods. The former two methods require an 
explicit estimate of the signal subspace via an eigen-decomposition of 
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the input covariance matrix which is extremely complex. Although the 
latter technique is much simpler, it does not achieve near full rank 
performance when the filter rank D is significantly less than the full rank 
N. 

5 Other known reduced rank techniques include a multi-stage 

Wiener filter as proposed by Goldstein et al. in "A Multistage 
Representation Of The Wiener Filter Based On Orthogonal Projections" 
IEEE Trans. Inform. Theory, 44 (7), November, 1998 and an adaptive 
interference suppression algorithm as proposed by Honig et al. in 
10 "Adaptive Reduced-Rank Residual Correlation Algorithms For DS- 

CDMA Interference Suppression" In Proc. 32 Asilomar Conf. Signals, 
Syst. Comput., November, 1998. These methods perform well and do 
not require an eigen-decomposition. The present invention is an 
improvement of these latter two methods. 

15 

BRIEF SUMMARY OF THE INVENTION 

In accordance with the present invention, the disadvantages of 
prior adaptive digital filters noted above have been overcome. The 

20 reduced rank adaptive filter of the present invention projects a received 

signal onto a lower dimensional subspace where the reduced rank 
subspace is iteratively constructed by multiplying the last basis vector by 
the received sample covariance matrix. The reduced rank filter of the 
present invention achieves full rank performance with a projected signal 

25 subspace having a significantly smaller dimension than the dimension of 

the received signal subspace. Moreover, for interference suppression 
applications, the optimum filter rank does not substantially increase with 
the dimension of the received signal subspace. As such, the method of 
the present invention enables faster tracking and convergence with 

30 significantly less training samples than can be achieved with prior 

techniques. 
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More particularly, the method of the present invention filters 
successive, received signal samples to provide a desired signal wherein a 
group of N successive samples form a received sample vector of digital 
data having a N x 1 dimension. The method includes generating D + 1 
5 basis vectors where D is less than N and represents the reduced rank 

dimension. Each successive basis vector is generated from a given or an 
estimated steering vector and successively greater powers of a 
covariance matrix for a sequence of received sample vectors of data. 
The first basis vector is formed from the steering vector. A reduced rank 

10 vector of digital data having a D x 1 dimension is generated from a 

matrix of D basis vectors and the received sample vector of data. D 
filter coefficients are generated from correlations between pairs of the 
basis vectors. The desired signal output from the filter is generated from 
the filter coefficients and the reduced rank vector of data. 

15 The method of the present invention can be used in any adaptive 

filtering application including echo and noise cancellation, channel 
equalization, radar processing, adaptive antenna arrays and interference 
suppression. The method of the present invention is particularly 
advantageous for applications where long filter lengths are required and 

20 fast convergence and tracking are important. One such application is for 

space-time interference suppression in wireless Code-Division Multiple 
Access (CDMA) systems. 

These and other advantages and novel features of the present 
invention, as well as details of an illustrated embodiment thereof, will be 

25 more fully understood from the following description and drawings. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS 
OF THE DRAWING 

30 Fig. 1 is a block diagram of the reduced rank adaptive filter of 

the present invention for a wireless CDMA communication system; 
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Fig. 2 is a flow chart illustrating the steps of the reduced rank 
adaptive digital filtering technique of Fig. 1; 

Fig. 3 is a block diagram illustrating the generation of the 
received signal sample covariance matrix and steering vector; 
5 Fig. 4 is a diagram illustrating the generation of the basis vectors; 

and 

Fig. 5 is a flow chart illustrating the method of selecting the 
dimension of the reduced rank filter. 

10 DETAILED DESCRIPTION OF THE INVENTION 

The reduced rank adaptive digital filtering technique 10 of the 
present invention is illustrated in Fig. 1 for interference suppression in a 
wireless CDMA communication system having a receiver 12. Although 

15 illustrated for a wireless communication system, the filtering technique 

of the present invention can be used in any adaptive digital filter 
application. The receiver 12 picks up a received signal via an antenna 
14 coupled to a radio frequency (RF) amplifier and mixer 16. The 
output of the amplifier/mixer 16 is coupled to an analog to digital 

20 converter 18 to provide a serial stream of baseband digital data samples 

representing the received signal. The output of the digital to analog 
converter 18 is coupled to a serial to parallel converter 20. The serial to 
parallel converter 20 takes N serial samples and outputs the N samples 
in parallel to provide a received sample vector r(i) of digital data 

25 having a dimension of N x 1 . 

The received sample vector r(i), which is referred to as an 
observation signal, is coupled to an adaptive digital filter 22, the output 
of which is an estimate of the desired signal b(i) , referred to as the 

approximate desired signal b(i) . A generator 24 produces a matrix M of 
30 basis vectors with a NxD dimension where D is less than N and 

represents the reduced rank dimension. The generator 24 also produces 
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a vector of filter coefficients c of dimension D x 1 from a sequence 
received sample vectors r(z) for i = l,...,B and from a sequence of the 
decision derived desired signals b(f), i — 1, . . . , B or a sequence of given 
desired signals or training symbols b(i), i = 1, B . The given training 
5 symbols b(f) initially used by the generator 24 may be stored, for 

example, in a receiver memory (not shown). 

The adaptive filter 22 generates a reduced rank or projected 
received signal vector 7(i) of dimension D x 1 by multiplying the 
Hermitian transpose of the matrix M by the received sample vector r(i) 
10 as follows. 

r(i)=M*r(i) 

The adaptive filter 22 generates the approximate desired signal b(i) by 
multiplying the Hermitian transpose of the filter coefficient vector by the 
reduced rank vector as follows. 

15 b(i) = c i r(i) 

The approximate desired signal b(i) is applied to a conventional sheer 
23 that essentially rounds the approximate desired signal b(i) to a 
desired signal level b(i) which may be used to generate the filter 
coefficients in a decision directed training mode. 

20 The generator 24 generates D + 1 basis vectors where each 

successive basis vector is generated from a given or an estimated 
steering vector p and successively greater powers of a sample covariance 
matrix R for a sequence of B received sample vectors r(i), i = l,...,B. 
The first basis vector is formed from the steering vector p as described 

25 in detail below with reference to Fig. 4. The generator 24 generates a 

D x 1 vector of filter coefficients c from a matrix and vector of 
correlations between pairs of generated basis vectors. 
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The generator 24 and adaptive digital filter 22 operate in 
accordance with the method depicted in Fig. 2. At block 26, the 
generator 24 determines the sample covariance matrix R for the 
sequence of received sample vectors r(i), i = l,...,B and the steering 
5 vector p from the sequence of received sample vectors r(i), i = 1, B, 

and the sequence of training/desired signals b(i), i = l, ... B in 
accordance with the block diagram of Fig. 3. Specifically, at block 28, 
the Hermitian transpose of the received sample vector is formed and at 
30 is multiplied by the received sample vector and a weighting factor of 

10 1 - w. The resulting product is added at 32 to the product formed at 34 

by multiplying the weighting factor w by the previous covariance matrix 
so as to form the sample covariance matrix R t for the sequence of 
received sample vectors r(l), r(i) . At 36, the conjugate of the given 
or desired signal b(i) is formed and at 38 is multiplied by the received 

15 sample vector r(z) and the weighting factor 1 - w. The resulting 

product is added at 40 to the product of the previous steering vector p.^ 
and the weighting factor w to form the new steering vector p t . The 
steering vector p t generated from the sequence of given training 
symbols b(i), i = l,...,B or from the sequence of desired filter outputs 

20 b(i),i = l,...,B as shown in Fig. 3 is an example of an "estimated" 

steering vector. It should be appreciated that Fig. 3 is only one example 
of a technique for generating the estimated steering vector. Other 
known techniques for generating an estimated steering vector can be 
used as well. A "given" steering vector can be used in accordance with 

25 the present invention as well. A "given" steering vector is a steering 

vector that is already known to the receiver and may be stored, for 
example, in a memory or buffer of the receiver. 
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At block 44, the generator 24 generates D + 1 basis vectors v 0 
through v D . The basis vectors are generated in accordance with the 
diagram depicted in Fig. 4. The first or initial basis vector v 0 is set 
equal to the steering vector p. Each successive basis vector is formed by 
multiplying the immediately preceding basis vector by the covariance 
matrix R so that the basis vectors are generated as follows. 

v 0 =p, v i =Rp, v 2 = R 2 p,...,v DA = R D ~ 1 p, v D =R D p 
The generator 24 at block 46 forms the matrix M of basis vectors from 
the basis vectors v 0 through v D _ x so that the matrix M has a NxD 
dimension. Thereafter, the generator 24 determines the correlation 
scalars y x through y 2D _ x from the basis vectors v 0 through v D 
determined at block 44. The correlation scalar y l+J is generated by 
multiplying the Hermitian transpose of the basis vector v t by the basis 
vector Vj as follows. 

Yi=v 0 f v x , Yi=v Q \ 2 ...y D =v Q \ D , y D+l =v D v x , y D+2 =v D v 2 ...y 2D _ 1 =v D v 

Also y 0 =v o \ 0 and v 0 is assumed to be normalized to that y 0 =l. 

The generator 24 at block 54 forms a correlation matrix T from 
the correlation scalars y x through y 2D _ x as follows. 

>i Ti-'-Yd 

Yd Yd+i-'-Yid-i. 

At block 54 the generator 24 also generates a reduced rank correlation 
vector y which is a D x 1 vector formed of the correlation scalars y 0 
through y D _ x as follows. 
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To 




Jd-i. 

The generator 24 at block 56 then generates the filter coefficient vector 
c which is a D x 1 vector generated by solving the following set of 
linear equations. 

Tc = y 

The matrix M of basis vectors formed at block 46 and the filter 
coefficients c generated at block 56 are applied to the adaptive digital 
filter 22 so as to generate the reduced rank vector r(i) and the 
approximate desired signal b{i) output from the filter 22 in accordance 
with blocks 58 and 60. In particular, at block 58, the adaptive digital 
filter 22 generates the reduced rank vector r(i) by multiplying the 
Hermitian transpose of the matrix M of basis vectors by the received 
sample vector of digital data r(i) as follows. 

r(i)=MV(i) 

At block 60, the filter 22 generates the approximate desired signal b(i) 
by multiplying the Hermitian transpose of the D x 1 filter coefficient 
vector c by the reduced rank vector r (i) as follows. 

It has been found that the reduced rank adaptive digital filter 
described above achieves near full rank performance with a value of D 
equal to 8 or less for a Direct Sequence (DS) Code Division Multiple 
Access (CDMA) communication system. For other applications, the 
value of D can be determined in accordance with the method depicted in 
Fig. 5. As shown in Fig. 5, at block 60 an index n is initialized to 1. 
Thereafter, at block 62 the n th basis vector v n is generated by 
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multiplying the previous basis vector v n _ x by the covariance matrix R 
where the first basis vector is the steering vector p so that 
v„ = R"p 

as discussed above. Thereafter, a matrix M n _ x of basis vectors is formed 
at block 64 as follows. 

M n-l = k»V 1 ,...V n _ 1 ] 

At block 66, the orthogonal projection u n , which is a N x 1 vector, is 
generated as follows. 

u n = v n -M n _ 1 (M t n -iM n _ i y 1 M n _ 1 \ n 
The respective lengths 1 and g of the vectors u n and v n are computed at 
block 67 in accordance with the following equations. 

' = IKI 
* - IKI 

where 1 and g are scalars. Thereafter, at block 68, it is determined 
whether the length 1 divided by the length g is greater than a small 
constant 8 , where 5 may be set equal to 0.01 for example. If so, n is 
incremented by 1 and steps 62 through 68 are repeated. Steps 62 
through 70 are repeated until 1 divided by g is not greater than 5 for a 
given value of n. At that point, D is set equal to n - 1 . 

In an alternative method, the value of D is chosen to minimize 
the "a posteriori Least Squares cost function C" where 

C = Y\b{m)-c\r{mf 

m=m 0 

where c D is the reduced rank filter coefficient vector generated from the 
sequence of b(Y}, ...b(B) and the sequence of r(l), ... r(B) and B<m c . 

The reduced rank adaptive filter of the present invention projects 
the N x 1 received signal vector r(i) onto the lower dimensional 
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subspace which results in a D x 1 vector r (i) where D can be much less 
than N while still allowing the filter 10 to achieve near full rank 
performance. For interference suppression applications, the optimum 
filter rank D does not substantially increase with the dimension of the 

5 signal subspace N. This enables much faster tracking and convergence 

than can be achieved with prior techniques. Furthermore, the method of 
the present invention does not rely on an estimate of the signal space via, 
for example, eigen-decomposition so that the complexity of the present 
method is substantially reduced. Although the method was described in 

10 detail for a wireless communication system application, it can be used in 

any adaptive filtering application including, echo and noise cancellation, 
channel equalization, radar processing, adaptive antennae arrays, 
interference suppression, etc. 

Many modifications and variations of the present invention are 

15 possible in light of the above teachings. Thus, it is to be understood 

that, within the scope of the appended claims, the invention may be 
practiced otherwise than as described hereinabove. 

What is claimed and desired to be secured by Letters Patent is: 

20 
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CLAIMS 

1. A method of filtering successive received signal samples, 
a group of N successive samples forming a received sample vector of 

5 digital data having a N x 1 dimension to provide an approximate desired 

signal comprising: 

generating a set of basis vectors where each successive basis 
vector is a function of a given or an estimated steering vector and 
successively greater powers of a covariance matrix for a sequence of 
10 received sample vectors of data with the initial basis vector being 

formed from the steering vector; 

generating a reduced rank vector of digital data having a D x 1 
dimension, where D is less than N, from a matrix of D basis vectors and 
a received sample vector of data; 
15 generating a D x 1 filter coefficient vector from the generated 

basis vectors; and 

generating the approximate desired signal from the filter 
coefficients and the reduced rank vector of data. 

2. A method of filtering successive received signal samples 
as recited in claim 1 where D + 1 successive basis vectors are generated 
by multiplying an immediately preceding basis vector of data by the 
covariance matrix for the received sample vector. 

3. A method of filtering successive received signal samples 
as recited in claim 1 wherein the reduced rank vector is generated by 
multiplying the Hermitian transpose of the matrix of basis vectors by the 
received sample vector of data. 

4. A method of filtering successive received signal samples 
as recited in claim 1 including forming a D x 1 correlation vector with 
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correlation scalars, each correlation scaler y i formed by multiplying the 

Hermitian transpose of the initial basis vector by the i th basis vector for 
i equal 0 through D-l. 

5. A method of filtering successive received signal samples 
as recited in claim 4 including forming a D x D correlation matrix T 
from 2D-1 correlation scalars generated from pairs of the basis vectors 
and generating a D x 1 vector of filter coefficients c by solving a set of 
linear equations defined by Tc = y where y is the correlation vector. 

6. A method of filtering successive received signal samples 
as recited in claim 1 including generating the desired signal by 
multiplying a Hermitian transpose of the filter coefficient matrix by the 
reduced rank vector of data. 

7. A method of filtering successive received signal samples, 
a group of N successive samples forming a received sample vector of 
digital data having a N x 1 dimension to provide an approximate desired 
signal comprising: 

determining a received sample covariance matrix from a 
sequence of received sample vectors of digital data; 
determining a steering vector of data; 

generating a reduced rank vector of digital data having a D x 1 
dimension where D is less than N from the received sample vector of 
data and a matrix of D basis vectors of data wherein each successive 
basis vector of data is generated by multiplying an immediately 
preceding basis vector of data by the covariance matrix and the initial 
basis vector is formed from the steering vector of data; 

generating D filter coefficients from a plurality of correlations 
between pairs of basis vectors of data; and 
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generating the approximate desired signal from the filter 
coefficients and the reduced rank vector of data. 

8. A method of filtering successive received signal samples 
as recited in claim 7 wherein the reduced rank vector is generated by 
multiplying the Hermitian transpose of the matrix of basis vectors by the 
received sample vector of data. 

9. A method of filtering successive received signal samples 
as recited in claim 7 including forming a D x 1 correlation vector with 
correlation scalars, each correlation scaler y l formed by multiplying the 

Hermitian transpose of the initial basis vector by the i th basis vector for 
5 i equal 0 through D-l . 

10. A method of filtering successive received signal samples 
as recited in claim 9 including forming a D x D correlation matrix T 
from 2D-1 correlation scalars generated from pairs of the basis vectors 
and generating a D x 1 vector of filter coefficients c by solving a set of 

5 linear equations defined by Yc = y where y is the correlation vector. 

11. A method of filtering successive received signal samples 
as recited in claim 7 including generating the desired signal by 
multiplying a Hermitian transpose of the filter coefficient matrix by the 
reduced rank vector of data. 

12. A method of filtering successive received signal samples, 
a group of N successive samples forming a received sample vector of 
digital data having a N x 1 dimension to provide an approximate desired 
signal comprising: 
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generating D + 1 basis vectors wherein D is less than N and each 
successive basis vector is a function of a given or an estimated steering 
vector and successively greater powers of a covariance matrix for a 
sequence of the received sample vectors of data with the first basis 
vector being formed from the steering vector; 

generating a reduced rank vector of digital data having a D x 1 
dimension from a matrix of D basis vectors and a received sample vector 
of data; 

generating 2 D - 1 correlations between pairs of basis vectors; 

generating D filter coefficients from a matrix of said correlations 
and a vector of D correlations; and 

generating the approximate desired signal from the filter 
coefficients and the reduced rank vector of data. 

13. A method of filtering successive received signal samples 
as recited in claim 12 where D + 1 successive basis vectors are 
generated by multiplying an immediately preceding basis vector of data 
by the covariance matrix for the received sample vector. 

14. A method of filtering successive received signal samples 
as recited in claim 12 wherein the reduced rank vector is generated by 
multiplying the Hermitian transpose of the matrix of basis vectors by the 
received sample vector of data. 

15. A method of filtering successive received signal samples 
as recited in claim 12 including forming a D x 1 correlation vector with 
correlation scalars, each correlation scaler y i formed by multiplying the 
Hermitian transpose of the initial basis vector by the i th basis vector for 
i equal 0 through D-l. 
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16. A method of filtering successive received signal samples 
as recited in claim 15 including forming a D x D correlation matrix T 
from 2D-1 correlation scalars generated from pairs of the basis vectors 
and generating a D x 1 vector of filter coefficients c by solving a set of 
linear equations defined by Tc = y where y is the correlation vector. 

17. A method of filtering successive received signal samples 
as recited in claim 12 including generating the desired signal by 
multiplying a Hermitian transpose of the filter coefficient matrix by the 
reduced rank vector of data. 

18. A method of filtering successive received signal samples, 
a group of N successive samples forming a received sample vector of 
digital data having a N x 1 dimension to provide an approximate desired 
signal comprising: 

selecting a dimension D for a reduced rank vector onto which a 
received sample vector of data is projected, where D is less than N; 

generating a set of basis vectors wherein each successive basis 
vector is a function of a given or an estimated steering vector and 
successively greater powers of a covariance matrix a sequence of 
received sample vectors of data with the first basis vector being formed 
from the steering vector; 

generating a reduced rank vector of digital data having a D x 1 
dimension from a matrix of D basis vectors and the received sample 
vector of data; 

generating D filter coefficients from the generated basis vectors; 

and 

generating the approximate desired signal from the filter 
coefficients and the reduced rank vector of data. 
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19. A method of filtering successive received signal samples 
as recited in claim 18 where D is selected to be less than or equal to 8. 

20. A method of filtering successive received signal samples 
as recited in claim 18 where D is selected to minimize an a posteriori 
Least Squares cost function. 

21. A method of filtering successive received signal samples 
as recited in claim 20 where D + 1 successive basis vectors are 
generated by multiplying an immediately preceding basis vector of data 
by the covariance matrix for the received sample vector. 

22. A method of filtering successive received signal samples 
as recited in claim 20 wherein the reduced rank vector is generated by 
multiplying the Hermitian transpose of the matrix of basis vectors by the 
received sample vector of data. 

23. A method of filtering successive received signal samples 
as recited in claim 20 including forming a D x 1 correlation vector with 
correlation scalars, each correlation scaler y t formed by multiplying the 
Hermitian transpose of the initial basis vector by the i th basis vector for 
i equal 0 through D-l. 

24. A method of filtering successive received signal samples 
as recited in claim 23 including forming a D x D correlation matrix T 
from 2D-1 correlation scalars generated from pairs of the basis vectors 
and generating a D x 1 vector of filter coefficients c by solving a set of 
linear equations defined by Tc = y where y is the correlation vector. 
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25. A method of filtering successive received signal samples 
as recited in claim 20 including generating the desired signal by 
multiplying a Hermitian transpose of the filter coefficient matrix by the 
reduced rank vector of data. 
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ABSTRACT OF THE DISCLOSURE 

A reduced rank adaptive digital filtering method is described for 
a received signal consisting of a sequence of N x 1 received vectors. 
Each received vector is formed from a group of N successive samples. 
D + 1 basis vectors are generated where D is less than N and the 
dimension of a desired reduced rank subspace. Each successive basis 
vector is generated by multiplying an immediate preceding basis vector 
by the covariance matrix for the sequence of received sample vectors 
and the first basis vector is formed from a given or estimated steering 
vector. D filter coefficients are generated from correlations between 
pairs of basis vectors. The adaptive digital filter of the present invention 
achieves near optimal rank performance with substantially fewer 
training symbols than heretofore possible. 
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